Dynamics of alliance formation and the egalitarian revolution 



Sergey Gavrilets*^ Edgar A. Duenez-Guzman^ and Michael D. Vose^ 
Departments of *Ecology and Evolutionary Biology, ^Mathematics, 
and ^Computer Science, University of Tennessee, Knoxville, TN 37996, USA. 
*corresponding author. Phone: 865-974-8136, fax: 865-974-3067, 
email: sergey@tiem.utk.edu 

September 1, 2008 
Abstract 

Background. Arguably the most influential force in human history is the formation of social coalitions and 
alliances (i.e., long-lasting coalitions) and their impact on individual power. Understanding the dynamics of 
alliance formation and its consequences for biological, social, and cultural evolution is a formidable theoretical 
challenge. In most great ape species, coalitions occur at individual and group levels and among both kin and non- 
kin. Nonetheless, ape societies remain essentially hierarchical, and coalitions rarely weaken social inequality. 
In contrast, human hunter-gatherers show a remarkable tendency to egalitarianism, and human coalitions and 
alliances occur not only among individuals and groups, but also among groups of groups. These observations 
suggest that the evolutionary dynamics of human coalitions can only be understood in the context of social 
networks and cognitive evolution. 

Methodology/Principal Findings. Here, we develop a stochastic model describing the emergence of net- 
works of allies resulting from within-group competition for status or mates between individuals utilizing dyadic 
information. The model shows that alliances often emerge in a phase transition-like fashion if the group size, 
awareness, aggressiveness, and persuasiveness of individuals are large and the decay rate of individual affinities 
is small. With cultural inheritance of social networks, a single leveling alliance including all group members can 
emerge in several generations. 

Conclusions/Significance. We propose a simple and flexible theoretical approach for studying the dynamics 
of alliance emergence applicable where game-theoretic methods are not practical. Our approach is both scalable 
and expandable. It is scalable in that it can be generalized to larger groups, or groups of groups. It is expandable 
in that it allows for inclusion of additional factors such as behavioral, genetic, social, and cultural features. 
Our results suggest that a rapid transition from a hierarchical society of great apes to an egalitarian society of 
hunter-gatherers (often referred to as "egalitarian revolution") could indeed follow an increase in human cognitive 



1 



abilities. The establishment of stable group-wide egalitarian alliances creates conditions promoting the origin of 
cultural norms favoring the group interests over those of individuals. 
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Introduction 

Coalitions and alliances (i.e., long-lasting coalitions) are often observed in a number of mammals including hye- 

fl). 



nas, wolves, lions, cheetahs, coatis, meerkats, and dolphins 



In primates, both kin and non-kin, and both 



within-group and group-level coalitions are a very powerful means of achieving increased reproductive success 

"lar" 



via increased dominance status and access to mates and other resources (Ul|2l,|3l|J,|5tlaLZ|). In humans, coalitions 
occurs at many different levels (ranging from within-family to between-nation states) and represent probably the 
most dominant factor in social interactions that has shaped human history 

The evolutionary forces emerging from coalitionary interactions may have been extremely important for the 
origin of our "uniquely unique" species HQ). For example, it has been argued that the evolution of human brain 
size and intelligence during Pleistocene was largely driven by selective forces arising from intense competition 
between individuals for increased social and reproductive success (the "social brain" hypothesis, also known as the 



and reproductive success (the social bra 
s; 
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271)). Coalition formation is 



"Machiavellian intelligence" hypothesis; 
one of the most powerful strategies in competitive interactions and thus it should have been an important ingredient 
of selective forces acting in early humans. Moreover, one can view language as a tool that originally emerged for 
simplifying the formation and improving the efficiency of coalitions and alliances. It has also been argued that 
the estabUshment of stable group- wide egalitarian alliances in early human groups should have created conditions 
promoting the origin of conscience, moralistic aggression, altruism, and other norms favoring the group interests 
over those of individual s (l28h. Increasing within-group cohesion should also promote the group efficiency in 
between-group conflicts 1291 Isoi and intensify cultural group selection Isii . 

In spite of their importance for biological, social and cultural evolution, our understanding of how coalitions 
and alliances are formed, maintained and break down is limited. Existing theoretical approaches for s tudying 



coalitions in animals are deeply rooted in cooperative game theory, economics, and operations research (132 , 
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Is^lssl). These approaches are usually limited by consideration of coalitions of two individuals against one, focus 
on conditions under which certain coalitions are successful and/or profitable and assume (implicitly or explici tly) 
that individuals are able to evaluate these conditions and join freely coalitions that maximize their success ([36 , 

ally do not capture the dynamic nature of coalitions and/or are 
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not directly applicable to individuals lacking the abilities to enter into binding agreements and to obtain, process, 
and use complex information on costs, benefits, and consequences of different actions involving multiple parties 
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(Q). These approaches do not account for the effects of friendship and the memory of past events and acts which 
all are important in coalition formation and maintenance. Other studies emphasize the importance of Prisoner's 
Dilemma as a paradigm for the emergence of cooperative behavior in groups engaged in the public goods game 
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471) . These studies have been highly successful in identifying conditions that favor the evolution of cooperation 
among unrelated individuals in the face of incentives to cheat. Prisoner's Dilemma however is often not appropriate 
for studying coaUtionary behavior (1481 l49h especially when individuals cooperate to compete directly with other 
individuals or coalitions and within-coalition interactions are mutualistic rather than altruistic and the 

benefit of cooperation is immediate. The social network dynamics that result from coalition formation remain 
largely unexplored. 

Here, we propose a simple and flexible theoretical approach for studying the dynamics of alliance emergence 
applicable where game-theoretic methods are not practical. Our method is related to recent models of social 
network formation and games on graphs with dynamic linking 

(0 Q, H H Q H). In our novel approach, 
alliances are defined in a natural way (via affinity matrices; see below) and emerge from low-level processes. 



The approach is both scalable and expandable. It is scalable in that it can be generahzed to larger 



groups, or 



groups of groups, and potentially applied to modeling the origin and evolution of states 

571) . It is expandable in that it allows for inclusion of additional factors such as behavioral, genetic, social, and 



cultural features. One particular application of our approach is an analysis of conditions under which intense 
competition for a limiting resource between individuals with intrinsically different fighting abilities could lead to 
the emergence of a single levehng alliance including all members of the group. This application is relevant with 
regard to recent discussions of "egalitarian revolution" (i.e. a rapid transition from a hierarchical society of great 
apes to an egalitarian society of human hunter-gatherers, and whether it could have been triggered by an 

increase in human cognitive abihties 



Model 

We consider a group of N individuals continuously engaged in competition for status and/or access to a limited 
resource. Individuals differ with regard to their fighting abilities s; (1 < i < N). The attitude of individual 
i to individual j is described by a variable Xij which we call affinity. We allow for both positive and negative 
affinities. Individual affinities control the probabilities of getting coalitionary support (see below). The group state 
is characterized by an x iV matrix with elements which we will call the affinity matrix. 

Time is continuous. Below we say that an event occurs at rate r if the probability of this event during a short 
time interval dt is rdt. 

We assume that each individual gets engaged in a conflict with another randomly chosen individual at rate a 
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which we treat as a constant for simplicity. Each other member of the group is aware of the conflict with a constant 
probabiHty lo. Each individual, say individual k, aware of a conflict between individuals i and j ("initiators"), 
evaluates a randomly chosen initiator of a conflict, say, individual i, and helps him or not with probabilities h^i 
and 1 — hki, respectively. In the latter case, individual k then evaluates the other initiator of the conflict and helps 
him or not with probabilities h^j and 1 — h^j, respectively. We note that the coahtionary support may be vocal 
rather than physical dSSh . Below we will graphically illustrate the group state using matrices with elements hij 
which we will call interference matrices. 

The interference probabilities hij are given by an S-shaped function of affinity xij and are scaled by two pa- 
rameters: (3 and rj. A baseline interference rate j3 controls the probability of interference on behalf of an individual 
the affinity towards whom is zero; j3 can be viewed as a measure of individual aggressiveness (i.e., the readiness to 
interfere in a conflict) or persuasiveness (i.e., the ability to attract help). A slope parameter 77 controls how rapidly 
the probability of interference increases with affinity. In numerical simulations we will use function 



hk 



ki 



1-13 , 

1 H ^ exp(-77a;fei 

P 



Note that the probability of help hij changes from to (3 to 1 as affinity Xij changes from large negative values to 
zero to large positive values. 

For simplicity, we assume that interference decisions are not affected by who else is interfering and on which 
side. We also assume that individuals join coalitions without regard to their probability of winning. This assump- 
tion is sensible as a first step because predicting the outcomes of conflicts involving multiple participants and 
changing alliances would be very challenging for apes and hunter-gatherers. 

As a result of interference, an initially dyadic conflict may transform into a conflict between two coalitions. 
[Here, coalition is a group of individuals on the same side of a particular conflict.] The fighting ability Sj of a 
coahtion / with n participants is defined as s„n^, where s„ is the avera ge fi ghting ability of the participants. This 



6lh which captures a larger importance of 



formulation follows the classical Lanchester-Osipov square law js^Q, j 
the size of the coalition over the individual strengths of its participants. The probability that coalition / prevails 
over coalition J is Si/ {Si + Sj). 

Following a conflict resolution we update the affinities of all parties involved by a process analogous to rein- 
forcement learning ( l62h . The affinities of winners are changed by 5ww, of the losers by Su, the affinities of winners 
to losers by S^i, and those of losers to winners by (5;^. The (5-values reflect the effects of the costs and benefits 
of interference on future actions. It is natural to assume that the affinities of winners increase (Siu^ > 0) and 
those of antagonists decrease (S^i < 0, (J^u < 0). The change in the affinities of losers Su can be of either sign 
or zero. Parameters Sww, Swi, Siw and 6u are considered to be constant. We note that a negative impact of costs 
of interfering in a conflict on the probability of future interferences can be captured by additionally reducing all 
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Table 1: Main dinamic variables, parameters, and other variables, functions, and statistics. 



Xij 


affinity of individual i to individual j 


N 


group size 


Si 


fighting ability of individual i 


a 


conflict initiation rate 




awareness 


P 


baseline interference rate 


V 


slope parameter 




changes in affinity after conflict resolution 




affinity decay rate 


K 


strength of social network inheritance 


7 


birth rate 


hij 

Si = sn2 
Si /{Si + Sj) 
X, 

= Efc bk/A, 
Hx,Hy 


probability that individual i helps individual j; is given by 
an S- shaped function of affinity with parameters (3 and rj 
strength of coalition / with n members and average fighting ability s 
probability that coalition / wins a conflict with coalition J 
proportion of conflicts won by individual i since birth 
expected social success of individual i; Ai is the age 
of individual i and hk is the benefit of the fcth conflict 
standard deviations of Xi and Yi in the group (measures of inequality) 
clustering coefficients and the average probability of help in an alliance 



affinities between members of a coalition by a fixed value 5. 

We assume that coalitions are formed and conflicts are resolved on a time-scale much faster than that of conflict 
initiation. Finally, to reflect a reduced importance of past events relative to more recent events in controlling one's 
affinities, affinities decay towards at a constant rate ji (l63h . Table 1 summarizes our notaion. 



Results and Their Biological Interpretation 

To gain intuition about the model's behavior we ran numerical simulations with all affinities initially zero. We 
analyzed the structure of the interference matrix hij, looking for emerging alliances. We say individuals i and j 
are allies if their interference probabilities hij and hji both exceed the baseline interference rate (3 by at least 50%. 
An alliance is a connected network of allies. 
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We also measured a number of statistics including the average and variance of affinities, the proportion of 
individuals who belong to an alliance, the number and sizes of alhances, the clustering coefficients C*^^^ and C'^^^ 



(1641) . related to the probability that two allies of an individual are themselves allies. The average interference 
probability and the clustering coefficients can be interpreted as measuring the "strength" of alliances. 

To make interpretation of model dynamics easier, we computed the proportion Xi of conflicts won since birth, 
and the expected social success Yi = bk/Ai, where Ai is the age of individual i, the sum is over all conflicts 
k he has participated in, the benefit 6^ is 1/ Uk if i was a member of a winning coalition of Uk individuals, and 
bk is if i was on the losing side. Although in our model the probability of winning always increases with the 
coalition size, the benefit bk always decreases with the coalition size. The net effect of the alliance size on the 
expected benefits of its members will depend on the sizes and composition of all alliances in the group. Note that 
our interpretation of Yi as a measure of expected social success makes sense both if all members on the winning 
side share equally the reward or if the spoils of each particular conflict goes to a randomly chosen member of the 
winning coalition. The former may be the case when the reward is an increase in status or rank. The latter may 
correspond to situations similar to those in baboons fighting over females, where members of the winning coahtion 



may race to the female and whoever reaches her first becomes the undisputed consort for some time (14 8h . Nonequal 
sharing of benefits can be incorporated in the model in a straightforward way. Note also that being a member of a 
losing coalition always reduces relative social success. 

We also calculated the standard deviations Hx and Hy of Xi and Yi values. These statistics measure the 
degree of "social inequahty" in the group. 

Figure [U illustrates some coalitionary regimes observed in simulations using a default set of parameters (a = 
l,f3 = 0.05, Syjyj — 1, Sii = 0.5, Syji = " 0.5, — —0.5, r] = 0.5, oj = 0.5, fia = 0.05) unless noted otherwise. 
This figure shows the N x N interference matrices using small squares arranged in an x array with each 
of the squares color-coding for the corresponding value of hij using the gray scale. The squares on the diagonal 
are painted black for convenience. In all examples, individual strengths Si are chosen randomly and independently 
from a uniform distribution on [0, 10] resulting in strong between-individual variation. 

Emergence of alliances. In our model, the affinity between any two individuals is reinforced if they are on 
a winning side of a conflict and is decreased if they are on the opposite sides; all affinities also decay to zero 
at a constant rate. The resulting state represents a balance between factors increasing and decreasing affinities. 
Although the emergence of alliances is in no way automatic, simulations show that under certain conditions they 
do emerge. The size, strength, and temporal stability of alhances depend on parameters and may vary dramatically 
from one run to another even with the same parameters. However, once one or more alliances with high values 
of C'(i),C(2) and ft: are formed, they are typically stable. Individuals belonging to the same alliance have very 
similar social success which is only weakly correlated with their fighting abilities. That is, the social success is 
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(a) 



(c) 



(e) 



(b) 




(d) 




(f) 



Figure 1: Interference matrices at time 1000. Values of hij are gray-scale coded from (white) to 1 (black), with diagonal 
elements set to black. The smallest squares on the diagonal represent unaffiliated individuals. For display purposes, alliances are 
ordered according to their clustering coefficients C'^' so that stronger alliances occur first along the diagonal. Parameters have 
default values except where noted, (a) iV = 10. (b) iV = 20. (c) iV = 30. (d) TV = 20, 5u = -0.5. (e) N = 20, fi = 0.1. 
(f) iV = 20, tJ = 0.25. 



now defined not by the individual?s fighting ability but by the size and strength of the alliance he belongs to. 
Individuals from different alliances can have vastly different social success, so that the formation of coalitions and 
alhances does not necessarily reduce social inequality in the group as a whole. 

Phase transition. We performed a detailed numerical study of the effects of individual parameters of the prop- 
erties of the system. As expected, increasing the frequency of interactions (which can be achieved by increasing 
the group size N, the awareness probability uj, baseline interference rate /5, or the slope parameter 77) and reducing 
the affinity decay rate fi all promote alliance formation. Most interestingly, some characteristics change in a phase 
transition-like pattern as some parameters undergo small changes. For example. Figure |2] show that increasing 
N, u!, f3, 77, or decreasing n result in a sudden transition from no alliances to at least one very strong alliance with 
all members always supporting each other. Parameter 6u has a similar but less extreme effect, whereas parameters 
Siui and 6iw have relatively weak effects (Supplementary Information). Similar threshold-hke behavior is exhibited 
by the C^^) -measure, the average probabihty of help h within the largest alliance, the number of alliances, and the 
numbers of alliances with C^^^ > 0.5 and with h > 0.5 (Supplementary Information). Interestingly, formation of 
multiple alliances is hindered when affinities between individuals fighting on the same side decrease as a result of 
losing (i.e., if Sn < 0). 

Cultural inheritance of social networks. Next, we extended the model to larger temporal scales by allowing 
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Figure 2: Tiikey plots for the effects of A*', a;, /3, rj, ^, Sn on the C'^' measure of the largest alliance. Each graph shows 
the effect of changing a single parameter from its default value (results for each parameter value are averaged over 20 runs, 
using data from time 1000 to 2000). The vertical lines extend from minimum to maximum observations, the dashed lines depict 
averages, and the boxes extend from lower to upper quartiles. 



for birth/death events, and the cultural inheritance of social networks. New individuals are born at a constant rate 7. 
Each buth causes the death of a different randomly chosen individual. We explored two rather different scenarios 
of cultural inheritance. In the first, the offspring inherits the social network of its parent who is chosen among all 
individuals with a probability proportional to the rate of social success Yi. This scenario requires special social 
bonds between parents and offspring. In the second, each new individual inherits affinities of its "role model" 
(chosen from the whole group either with a uniform probability or with a probability proportional to the rate of 
social success Yi). Under both scenarios, if individual i* is an offspring (biological in the first scenario or cultural 
in the second scenario) of individual i, then we set Xi*j — KXij for each other individual j in the group (parameter 
< K < 1 controls the strength of social network inheritance). In the parent-offspring case, the affinities of other 
individuals to the son are proportional to those to the father: Xji<- ~ nxji and Xi*i — xu* is set to k times the 
maximum existing affinity in the group. In the role model case, other individuals initially have zero affinities to 
the new member of the group: Xji- — 0. 

Stochastic equilibrium. If cultural inheritance of social networks is weak {k is small), a small number of 



8 



t= 100Q t = 200Q t = 3000 1 = 4000 t = 5000 t = 6000 t = 7000 t = 8000 
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t = 9000 t = 10000 




1 = 11000 t = 120D0 t = 130D0 t = 140D0 1 = 15000 1 = 16000 1 = 17000 



1 = 19000 1 = 20000 



Figure 3: Dynamics of interference matrix; no cultural inheritance (k. = 0). Parameters values are default with A'^ = 20 
and 7 — 0.001 (average life span 1000). See the legend of Fig.l. 



alliances are maintained across generations in stochastic equilibrium (see Figure |3]l. This happens because the 
death of individuals tends to decrease the size of existing alliances while new individuals are initially unaffiliated 
and may form new affinities. This regime is similar to coalitionary structures recently identified in a community of 



wild chimpanzees in Uganda (0) and in populations of bottlenose dolphins in coastal waters of Western Australia 



(l65h and eastern Scotland ( l66h . 

Egalitarian state. If cultural inheritance of social networks is faithful (k is large), the dynamics become 
dramatically different due to intense selection between different alliances. Now the turnover of individuals creates 
conditions for growth of alliances. Larger alliances increase in size as a result of their members winning more 
conflicts, achieving higher social success, and parenting (biologically or culturally) more offspring who themselves 
become members of the paternal alliance. As a result of this positive feedback loop (analogous to that of positive 
frequency-dependent selection), the system exhibits a strong tendency towards approaching a state in which all 
members of the group belong to the same alliance and have very similar social success in spite of strong variation 
in their fighting abilities. Figure 4 contrasts the egalitarian state with the stochastic equilibrium illustrated in Figure 
3 above. One can see that at the egalitarian state, the average affinity is increased while the standard deviation of 
affinity and the hierarchy measures are decreased. Although at the egalitarian state the correlation of individual 
strength and social success can be substantial, it does not result in social inequality. This "egalitarian" state can be 
reached in several generations. 

Cycling. However, the egalitarian state is not always stable. Under certain conditions the system continuously 
goes through cycles of increased and decreased cohesion (Figure [S^-c) in which the egalitarian state is gradually 
approached as one alliance eventually excludes all others. But once the egalitarian state is established (in Figure|5}l, 
around time 5200), it quickly disintegrates because of internal conflicts between members of the winning alliance. 
Figure |5}l illustrates one such cycle, showing that the dominant alliance remains relatively stable as long as the 
group excludes at least one member ("outsider"). 

Analytical approaches. Simple "mean-field" approximations help to understand model dynamics. These ap- 
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Figure 4: Contrasting a state with a small number of alliances maintained in stochastic equilibrium with an 

egalitarian state. The left column of graphs correspond to the run shown in Figure [3] which resulted in a small 

number of alliances maintained in stochastic equilibrium. The right column of graphs correspond to a run with 

K = 1 (complete cultural inheritance of social network) and fi — 0.025 (increased memory of past events) which 

resulted in an egalitarian regime. With several alliances present simultaneously (left graphs), the average affinity is 

small, the variance of affinities is large, the measures of social inequality and Hy are large, and the correlation 

between social success Yi and individual fighting ability Si is small. In the egalitarian state (right graphs), the 

average affinity is large, the variance of affinities is small, the measures of social inequality and Hy are small, 

and the correlation between social success Yi and individual fighting ability is large. 
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Figure 5: A coalitionary cycle; complete cultural inheritance (k = 1). Other parameters are as in Figure [5] (a) Average 
(blue) and standard deviation (red) of affinities, (b) Number of alliances (blue) and clustering coefficient for the largest 
alliance (red), (c) Proportions of individuals belonging to an alliance (red) and to the largest alliance (blue), (d) Dynamics of 
the interference matrix between time 4100 and 7200. 



proximations focus on the average a and variance v of affinities computed over particular coalitions (Supplemen- 
tary Information). For example, at an egalitarian state when all individuals have very high affinity to each other, the 
dynamics of a and v are predicted to evolve to particular stochastic equilibrium values, a* and v* . The egalitarian 
state is stable if the fluctuations of pairwise affinities around a* do not result in negative affinities. We conjecture 
that the egalitarian state is stable if a* > which is roughly equivalent to {a*)"^ > lOv*, which can be 

rewritten as 



2aNL0^ 



> 10 



var S 



+ 1-UJ' 



Here the mean 5 and variance var 5 are computed over the four ^-coefficients. Both the approximations and numer- 
ical simulations suggest that the egalitarian state cannot be stable with negative S. Increasing the population size 
N, awareness ui, average 6, and decreasing the affinity decay rate fi and variance var S all promote stability of the 
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egalitarian state. The agreement of numerical simulations with analytical approximations is very good given the 
stochastic nature of the process. Similar approximations can be developed for other regimes. In particular, one can 
show (Supplementary Information) that the stabihzing effect of "outsiders" on the persistence of alliances is espe- 
cially strong in small groups. This happens because successful conflicts against outsiders simultaneously increase 
the average a and decrease the variance v of the within-alliance affinities, with both effects being proportional to 
l/N. 



Discussion 

The overall goal of this paper was to develop a flexible theoretical framework for describing the emergence of 
aUiances of individuals able to overcome the power of alpha-types in a population and to study the dynamics and 
consequences of these processes. We considered a group of individuals competing for rank and/or some limiting 
resource (e.g., mates). We assumed that individuals varied strongly in their fighting abihties. If all conflicts 
were exclusively dyadic, a hierarchy would emerge with a few strongest individuals getting most of the resource 
However there is also a tendency (very smaU initially) for individuals to interfere in an ongoing 
dyadic conflict thus biasing its outcome one way or another Positive outcomes of such interferences increase the 
affinities between individuals while negative outcomes decrease them. Using a minimum set of assumptions about 
cognitive abilities of individuals, we looked for conditions under which long-lasting coalitions (i.e. alliances) 
emerge in the group. We showed that such an outcome is promoted by increasing the frequency of interactions 
(which can be achieved in a number of ways) and decreasing the affinity decay rate. Most interestingly, the model 
shows that the shift from a state with no alliances to one or more alliances typically occurs in a phase-transition 
like fashion. Even more surprisingly, under certain conditions (that include some cultural inheritance of social 
networks) a single alliance comprising all members of the group can emerge in which the resource is divided 
evenly. That is, the competition among nonequal individuals can paradoxically result in their eventual equality. 
We emphasize that in our model, egalitarianism emerges from political dynamics of intense competition between 
individuals for higher social and reproductive success rather than by environmental constraints, social structure, or 
cultural processes. In other words, within-group conflicts promote the buildup of a group-level alliance. In a sense, 
once alliances start to form, there is no other reasonable strategy but to join one, and once social networks become 
highly heritable, a single alliance including all group members is destined to emerge. 

Few clarifications are in order First, in our model coahtionary interactions are mutualistic in nature rather 
than altruistic. We note that there are not many examples of truly altruistic behavior outside of humans with 
some of those that were initially suggested to be altruistic under closer examination turning out to be kin-directed 
or mutualistic (HQ). 

Even in humans certain behaviors that are viewed as altruistic may have a rather different 
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origin. For example, food sharing may have originated as a way to avoid harassment, e.g. in the form of begging 
(45). In any case, modern human behavior is strongly shaped by evolved culture Isli and might not be a good 
indicator of factors acting during its origin. Second, in our model we avoided the crucial step of the dominant 
game-theoretic paradigm which is an explicit evaluation of costs and benefits of certain actions in controlling one's 
decisions. In our model, coalitions and alliances emerge from simple processes based on individuals using only 
limited "local" information (i.e., information on own affinities but not on other individuals' affinities) rather than 
as a solution to an optimization task. Our approach is justified not only by its mathematical simplicity but by 
biological realism as well. Indeed, solving the cost-benefit optimization tasks (which require rather sophisticated 
algebra in modern game-theoretic models) would be very difficult for apes and early humans especially given 
the multiplicity of behavioral choices and the dynamic nature of coalitions. Therefore treating coalitions and 
alliances in early human groups as an emergent property rather than an optimization task solution appears to 
be a much more realistic approach. We note that costs and benefits can be incorporated in our approach in a 
straightforward manner Third, one should be careful in applying our model to contemporary humans (whether 
members of modern societies or hunter-gathers). In contemporary humans, an individual's decision on joining 
coahtions will be strongly affected by his/her estimates of costs, benefits, and risks associated as well as by cultural 
beliefs and traditions. These are the factors explicitly left outside of our framework. 

Our results have implications for a number of questions related to human social evolution. The great apes' 
societies are very hierarchical; their social system is based on sharp status rivalry and depends on specific disposi- 
tions for dominance and submission. A major function of coahtions in apes is to maintain or change the dominance 
structure (jillsl); although leveling coalitions are sometimes observed (e.g., |^), they are typically of small size and 
short-lived. In sharp contrast, most known hunter-gatherer societies are egahtarian dsl l9lll0l). Their weak leaders 
merely assist a consensus-seeking process when the group needs to make decisions; at the band level, all main 
pohtical actors behave as equal. It has been argued that in egalitarian societies the pyramid of power is turned 
upside down with potential subordinates being able to express dominance because they find collective security in 
a large, group-wide political coalition (1^. One factor that may have promoted transition to an egahtarian soci- 
ety is the development of larger brains and better political/social intelligence in response to intense within-group 
competition for increased social and reproductive success Jl^flTilzsiuTl) . Our model supports these arguments. 
Indeed, increased cognitive abilities would allow humans to maintain larger group sizes, have higher awareness of 
ongoing conflicts, better abilities in attracting allies and building complex coalitions, and better memories of past 
events. The changes in each of these characteristics may have shifted the group across the phase boundary to the 
regime where the emergence of an egahtarian state becomes unavoidable. Similar effect would follow a change 
in mating system that would increase father-son social bonds, or an increase in fidelity of cultural inheritance of 
social networks. The fact that mother-daughter social bonds are often very strong suggests (everything else being 
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the same) females could more easily achieve egalitarian societies. The establishment of a stable group-wide egal- 
itarian alliance should create conditions promoting the origin of conscience, moralistic aggression, altruism, and 



other cultural norms favoring the group interests over those of individuals i2w . Increasing within-group cohesion 



will also promote the group efficiency in between-group conflicts (1291) and intensify cultural group selection. 

In humans, a secondary transition from egalitarian societies to hierarchical states took place as the first civ- 
ilizations were emerging. How can it be understood in terms of the model presented here? One can speculate 
that technological and cultural advances made the coalition size much less important in controlling the outcome 
of a conflict than the individuals' ability to directly control and use resources (e.g., weapons, information, food) 
that strongly influence conflict outcomes. In terms of our model, this would dramatically increase the variation 
in individual fighting abilities and simultaneously render the Lanchester-Osipov square law inapplicable, making 
egalitarianism unstable. 

Besides developing a novel and general approach for modeling coalitionary interactions and providing theoret- 
ical support to some controversial verbal arguments concerning social transitions during the origin of humans, the 
research presented here aUows one to make a number of testable predictions. In particular, our model has identified 
a number of factors (such as group size, the extent to which group members are aware of within-group conflicts, 
cognitive abilities, aggressiveness, persuasiveness, existence of outsiders, and the strength of parent-offspring so- 
cial bonds) which are predicted to increase the likelihood and size of alliances and affect in specific ways individual 
social success and the degree of within-group inequality. Existing data on coalitions in mammals (in particular, in 
dolphins and primates) and in human hunter-gatherer societies should be useful in testing these predictions and in 
refining our model. 
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Supporting Information 



Here, we present 

• some additional details on the computational methods used; 

• a set of figures (Figures [STTlSSI ) illustrating the effects of individual parameters on the coahtionary structure 
of the model achieved within a single generation; 

• a set of figures (Figures[S9land [ST0] i illustrating the effects of changes in multiple parameters simultaneously 
on the coalitionary structure of the model achieved within a single generation; 

• an outline of a mathematical method used to study the model analytically. 



Some details of computational methods 

Probabilities of help For an individual k aware of a conflict between individuals i and j, the probabilities of 
helping to i, to j, and of no interference are set to hf^i — hkihkj/2. h^j — hkihi.j/2 and (1 — hki){l ~ h^j), 
respectively. In numerical simulations, we set 



hki = 



IH —cxp[~r]Xkij 



where f3 and 77 are scaling parameters. Note that hki ~> 1 for Xki 00, hki ^ for Xki —00, and hki = P for 
Xki = 0. 

Numerical implementation The model dynamics were simulated using Gillespie's direct method (Gillespie 
1977). That is, the next event to happen is chosen according to the corresponding rates. The time interval until the 
next event is drawn from an exponential distribution with a parameter equal to the sum of the rates of all possible 
events. All rates are recomputed after each event. 



Reference 

• Gillespie, D. T. Exact stochastic simulation of coupled chemical reactions. Journal of Physical Chemistry 
81,2340-2361 (1977) 



Supplementary Figures and Legends 

Figures ISHISSI To obtain Fi gures ISTllS 8 1 we performed 20 runs for each parameter combination. Each of the 20 
runs was characterized by a single average value (computed over 100 observations taken between time 1000 to 
2000). All plots correspond to the Tukey Plots (i.e. show mean, min, max, quantile 1/4 and quantile 3/4), with 20 
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data points. Other parameters were set to default values (N = 20,5^^ = l,5u = 0.5, S^i = —0-5,Siw = — l,/? = 
0.05, /i = 0.1, 77 = 0.5,^^ = 0.5). 

Figures IS9l and lST0l To obtain Fi gure s [S9] and ISTOl we performed 40 runs for each parameter combination. 
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Figure S 1 : Effects of pai-ameters iV, (5;; , (5^,; , (5;^, , /3, /i, 77, w on the proportion of individuals belonging to a alliance 
for a default set of parameter values. 
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Figure S2: Effects of parameters N, 6u , 6wi,Siw, P, iJ;i],u} on the number of alliances for a default set of parameter 
values. 
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Figure S3: Effects of parameters N,Sii,Su]i,Siw, f3, fi,ri,uj on the size of the biggest alliance for a default set of 
parameter values. 
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Figure S4: Effects of parameters N, Su , S^i, Si^, P, ^,ri,uj on the C'^^^ measure of the largest alliance for a default 
set of parameter values. 
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Figure S5: Effects of parameters N, Sii,Swi,Siyj,P, ^, 77, w on the C*^^^ measure of the largest alliance for a default 
set of parameter values. 
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Figure S6: Effects of parameters N, Su , 5wi , 5iw , IJ-jIIjUj on the probability of help within of the largest alliance 
for a default set of parameter values. 
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Figure S7: Effects of parameters N, Su , 6^1 , Siw , P, fi, rj, u on the number of alliances with > 0.5 for a default 
set of parameter values. 
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Figure S8: Effects of parameters N, 6u, 5wi, Si^, (3, /i, rj, lo on the number of alliances with within-cluster proba- 
biUty of interference > 0.5 for a default set of parameter values. 
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Figure S9: Effects of parameters N, (3, /i, 77, uj on the number of individuals in alliances (first column) and the size 
of the largest alliance (second column) for (Ju,^, = 1.0, (5;; = 0.5,(5^; = —0.5,(5/^ = —1.0. First row: N = 10, 
second row: N = 20, third row: N = 30. 



30 



C^^^ of the largest alliance 



Number of alliances with C^^-* > 0.5 




N =10 





N = 30 

Figure SIO: Effects of parameters N, (3, fi, rj, uj on the C^^^ measure of the largest alliance (first column), and the 
number of alliances with C'^^ > 0.5 (second column) for (5^tu = 1.0, (5;; = 0.5, S^i = — 0.5, Si^ = —1.0. First 
row: N = 10, second row: = 20, third row: = 30. 
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Supplementary Methods: Mean field approximation for the dynamics of coalitions on the 
within-generation time-scale 

We consider a group of individuals in which conflicts occur at rate aN . Below we will use two types of 
averages: the average over a clique (i.e., a set of individuals who all are close allies), which we will denote as (^), 
and the average over all possible outcomes of the process, which we will denote as ^ or E{S,), where ^ is a random 
variable. 

Approximate dynamics of the mean and variance of affinities near an egalitarian state. We assume that 
all N individuals are close allies so that each individual aware of a conflict interferes in it. The average affinity of 
the group is 

After each conflict, each affinity value changes from Xij to Xij + where is a random variable describing the 
change in affinity of individual i to individual j. Let a = {x) be the expected average affinity. Since expectation 
and averaging are linear, the expected average affinity after a conflict can be written as 

a' = {x + e) = a + (e), 

All affinities continuously decay to at a constant rate /i. Therefore, the dynamics of a are described by a differ- 
ential equation 

^ = aNls) - aa. (SI) 
at 

Similarly, let v = (x^) — (a;)^ be the expectation of the variance taken over all possible outcomes of the 
process. Then the variance after a conflict is 

v' = i;(((x + e)2)-(x + e)2) 

= e[{x^) + 2{xe) + is') - {{xf + 2{x) {e) + (e)^)') 
= e({x') - {xf + is') - {sf) + 2{xe) - 2{x) (e) 
= v + W}-W- 

where, as an approximation, we assumed that x and e are independent with respect to the averaging operator, i.e., 

(xe) = (x) {e). 

All squares of affinities decay to at a constant rate 2ji. Therefore, the dynamics of v are described by a 
differential equation 

^=aN {{e^)-{£y)-2,,v. (S2) 
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First, we consider the expected change (e) in the affinity of a random pair of individuals after a conflict. There 
are three possibihties: 

• With probabiHty 1 / (^) , the two individuals are the initiators of the conflict. Since either of the two initiators 
can be on the winning side, the expected change in their affinity is 

f 5wL + 5lw 

On = . 

" 2 

Under our assumptions about the meaning of parameters, 5q is negative. 



• With probability 



2{N-2)/{ 



UJ, one of the two individuals is an "initiator" while the other was aware 
of the conflict and interfered on behalf of one side. Since there are four ways to distribute the two individuals 
over the winning and losing coalitions and each occurs with equal probability, the expected change in their 
affinity is 

-J Sww + SwL + Slw + Sll 
= : • 



• With probability 



l-l/(^)-2(iV-2)/(^) 



neither individual is the initiator of the conflict but both 



are aware of it and interfere in the conflict. The expected change in their affinity is 6. 
Therefore, 

(^) -2(iV-2) - 1 



1 , 2{N-2) ^ 
ivT oo H — w d - 



(^) 

2 



P5 



-{So-u;^ S). 



(S3a) 
(S3b) 



N{N- 1) ' ' N{N- 1) 
Then, equations ( ISllS3bl l predict that the average affinity in the egalitarian state evolves to an equilibrium value 



aN 



- , 4(7V-2) 



N{N- 1) 



w(l - uj)6 + 



N{N - 1) 



[So - S) 



(S4) 



The average affinity is positive only if S > 0. The last term in the brackets can be neglected relative to the first 
term even for small groups (e.g., N > 5). The second term in the brackets can be neglected for larger groups (e.g., 
N > 40) if Lu is not too small. Under these conditions, a* 
In a similar way and using the results above. 



^uj^ S. 



r2\ - 



2{N - 2) , ^{")- 2(1V - 2) 

UJ 02 -\ 



1 



(^) 



where 



r2 I r2 

Si- ^ , S2 



(^) 



A2 I X2 ,X2 , X2 

"ww ^ "WL + "Lty + "ll 



UJ^ 52, 



(S5) 
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More involved calculations show that 

w - — 



7V2(iV_ 1)2 ^ 



1 



1 4(jV - 2) (Ar-2)(Ar-3) 

7V(Af-l) ^"^iV(Af-l) ^"^ 7V(Ar-l) ^' 



where 



^1 

^2 
^3 



4Af(Ar-l)(Af-2) 
1 



N{N - l)(iV-2)(iV-3) 



(S6a) 
(S6b) 
(S6c) 
(S6d) 
(S6e) 

(S6f) 
(S6g) 
(S6h) 



»7^j k,l^i,j 

The term Ai can be interpreted as the expected value of A = SijEij + eijeji for a random pair of individuals 
(i and j). There are three cases to consider. 

• With probability 1 / (^) , the focal individuals are the initiators of the conflict. In this case, A = 26^. 



• With probability 



2(Ar-2)/Q 



CO, one of the two focal individuals is the initiator of the conflict while the 



other is aware of it. 



With probability 



1 - l/(^) - 2{N - 2)/(^) a;2^ both focal individuals are aware of the conflict. In the 



last two cases, A = Sq + 262 — 61. 
Therefore, 



N\ 



2(N-2) ,9 



2{N -2)- 1 



{6I + 262-61) 



(S7) 



(^) Q 

The term A2 can be interpreted as the expected value of A = eijEik + SijSkj + ^ij^jk + ^ij^ki for a random 
triple of individuals and k). There are three cases to consider. 



With probability 



3/( 



oj, two of the three focal individuals are the initiators of the conflict while the third 



is aware of it. In this case, A = {SSgS + 46 ) /3. 



34 



With probabiUty 3{N - 3)/ ('^) 
the two others are aware of it. 



w^, one of the three focal individuals is the initiator of the conflict while 



With probability 



l-3(7V-3)/(^)-3/Q 



(jj^, none of the three focal individuals are the initiators of 



the conflict but all are aware of it. 



To evaluate A in the last two cases, one needs to consider changes in affinities corresponding to all possible ways 
to assign three individuals to the winning and losing coalitions. This is done in the table below: 



winners 


losers 










ijk 




"WW 


"WW 


"WW 


"WW 


ij 


k 


SwW^WL 


SwwSlw 


SwwSwL 


SwwSlw 


ik 


3 


SwlSww 


S'WL 


SwlSlw 


SwlSww 


jk 


i 


Slw 


SlwSww 


SlwSww 


SlwSwl 


i 


jk 


S'WL 


6wl6ll 


SwlSll 


SwlSlw 


3 


ik 


SlwSll 


Slw 


SlwSwl 


SlwSll 


k 


ij 


SllSlw 


SllSwl 


SllSlw 


SllSwl 




ijk 


SIl 


SIl 


SIl 


SIl 



Using this table, 
1 



2<5 



wl + ^SwwSwL + ^SwwSlw + ^SllSlw + ^SllSwl + ^SlwSwl) 
=-s['^{Sww + Sll + 5lw + Swl)^ + '^■{Sww - SllY'] 

2 



-2 1 

=4(5 + -^{5ww - Sll) 



Therefore, 



4 

(T)- 



3(iV-3)^,2 , (^)-3(7V-3)-3 3' 



-co 



(45' + J [6ww - Sll^] . (S8) 



C.) (^) 

The term A3 can be interpreted as the expected value of A = SijEki for a random quartet of individuals ii,j, k 
and I). There are three cases to consider: 



• With probability 



6/(^) 



uj , two of the four focal individuals are the initiators of the conflict while the two 



others are aware of it. In this case, A = 63 = [45^ + 35o{6ww + Sll) + '^SllSww] /12 

With probability 4(iV - 4)/ (^) u)^, one of the four focal individuals is the initiator of the conflict while 



-^2 



the three others are aware of it. In this case, A = ^ . 



With probability 



l-4(iV-4)/(^)-6/(^) 



the conflict but aU are aware of it. In this case, A = ^ 



UJ , none of the three focal individuals are the initiators of 

:2 
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Therefore, 



^3 = -TTKUj'^h + 



(N\ 
, 2 , 



4(7V-4) 3 , (2)-4(iV-4)-6 4 

■W H TTJT^ W 



(S9) 



Keeping only the leading terms in {ej^) = {sij£ki) = (w^(5)^, which results in an equation for v: 



dv 
'di 



var^+ (1 -a;^)^ 



(SIO) 



where var S = S'^ — 6 . Higher order corrections (in 1/N) can be found in a straightforward way from the formula 
given above. 

Keeping only the leading terms in 1/N, the mean field approximation predicts the following equihbrium values 
at the egalitarian regime 

a = , 

/" 



YBX 5 + {I - UJ^)6 



2/x 

The egahtarian state is stable if the fluctuations of pairwise affinities around a* do not result in negative affini- 
ties. We conjecture that the egahtarian state is stable if a* > ^^/v*, which is roughly equivalent to (a*)^ > lOt;*, 
which in turn can be rewritten as 



> 10 



var 5 



+ 1-UJ' 



The strongest clique comprising A^^i individuals; other N2 = N — Ni individuals belong to weaker 
cliques. We assume that aU Ni individuals in the clique are close alhes that always help each other and never 
help outsiders. To evaluate the expected average over the chque (e), we need to find the expected value of A = Sij 
for a random pair from the strongest clique. One needs to consider five possibilities: 

• With probability 1 / (^) , the focal individuals are the initiators of the conflict. In this case, A = Sq. 

• With probability 2(7Vi - 2)a;/(^), one of the focal individuals is an initiator of a conflict involving another 
member of the chque while the other is aware of the conflict and interferes on behalf of one side. In this 
case, A = S. 

• With probability (^ 2^) / (?) ^^^^ both focal individuals are aware of and interfere in a conflict between 
two other members of the clique. In this case, A = 6. 

• With probability 2N2/{'^) w, one of the focal individuals is an initiator of a conflict involving an outsider 
while the other is aware of the conflict and interferes on behalf of the chque member. Assuming that the 
clique always wins, A = 6ww- 
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• With probability 



{Ni^2)N2l{ 



Lu , both focal individuals are aware of and interfere in a conflict be- 



tween a member of the clique and an outsider. Assuming that the clique always wins, A — 6ww- 
Therefore, 



2(7Vi - 2) 



ujS + 



2N- 



uj S + —j^LuSww + 



(iVi - 2)N, 



2. .2 



(^) 



UJ Sww ■ 



(Sll) 



(^) 

Assume that Ni ^ N ~ 1, N2 = I (i.e., the single outsider case). Then the dynamics of the average within-clique 
affinity a are described by equation 

(N-3\ 



da 



= aN 



So 



udww H 7/v\ — owiv 



.(^) (^) Q r.) Q 

Thus, the average affinity under the single outsider regime is predicted to evolve to 



lia. 



aN 



6{N~2) 

N{N- ly 



2{N ~3) 



N{N - 1) 



((5o + 2lj^,u,^) 



Keeping only terms of order OiXjN^ and larger in the brackets. 



aN 



- 4 - 2 



S) 



(S12) 



It is illuminating to compare this expression with expression (|S4| l approximating the average affinity under egali- 
tarian regime. Under the same assumptions, expression (IS4l i simplifies to 



aN 



LU^ 6 + — LU{1 - LU)S 



(SI 3) 



If N is not too large, a* can be substantially smaller than a* . It is in this situation when a single outsider can 
have a strong stabilizing effect on a small coalition. For example, let a = l,iV = 20, = 0.5,/^ = 0.05 and 
S^Yw = 1,Sll = 0.5, (Slvk = -0.5, Swl = -0.5 so that 6 = 0.125. Then a* = 15.00 but a* = 23.75, so that a 
single outsider significantly increases the average affinity of the clique. A single outsider will also reduce variance 
V, the effect of which will further strengthen the stability of the coalition. 
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